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Abstract. The crocodile lizard, Shinisaurus crocodilurus Ahl, 1930 , is a monotypic taxon, restricted in occurrence to south¬ 
ern China and northern Vietnam. Wild populations are presently suffering tremendous declines, mainly due to illegal 
poaching, habitat destruction, and fragmentation, which already led to the extinction of populations in Guangxi and 
Hunan provinces in China. In order to accelerate the discovery of so far unknown populations of S. crocodilurus and to 
identify suitable priority areas for conservation strategies, we determined the species’ potential distribution using correla¬ 
tive species distribution models (SDMs) based on locality records and a set of satellite-based environmental predictors. 
In addition, we evaluated the coverage of the species’ potential distribution with designated protected areas according to 
lUCN standards. The resulting SDM revealed potentially suitable habitats to be scattered and disconnected while being 
very small in size. Moreover, present coverage with nature reserves is extremely poor, underlining the urgent need for im¬ 
proved habitat protection measures and potential population restoration of S. crocodilurus. 

Keywords. Shinisauridae, Diploglossa, Conservation planning. Habitat suitability modelling. Species distribution model¬ 
ling, Southeast Asia, Vietnam. 


Introduction 

The crocodile lizard, Shinisaurus crocodilurus Ahl, 1930, 
is the only living representative of the monotypic family 
Shinisauridae, and despite its striking appearance, it was 
only described relatively recently (Hu et al. 1984, Zhang 
1991). The species usually is found along slow-flowing 
rocky streams in montane evergreen forests. The altitudi¬ 
nal range of this species was reported to reach from 200 
to 1,500 m in China and from 400 to 800 m in Vietnam. 
(Zhao et al. 1999, Le & Ziegler 2003, Huang et al. 2008). 
So far, the occurrence of S. crocodilurus has been confirmed 
from Guangxi and Guangdong provinces in southern Ghi- 
na while a couple of populations in Hunan and Guangxi 
provinces have already been extirpated (Huang et al. 2008, 
ZoLLWEG & Kuhne 2013, Z. Wu pers. comm.). In northern 
Vietnam, the species has been reported from the contigu¬ 


ous nature reserves Tay Yen Tu in Bac Giang Province and 
Yen Tu in Quang Ninh Province (Hecht et al., 2014, Le & 
Ziegler 2003, Ziegler et al. 2008, Nguyen et al. 2009). 
However, a variety of anthropogenic hazards have caused 
severe population declines within the last decades, reduc¬ 
ing estimated population densities in China from 6,000 to 
950 individuals between 1978 and 2008 (Zhao et al. 1999, 
Mo & Zou 2000, Huang et al. 2008). Illegal poaching for 
the international pet trade, traditional medicine and food 
represents the main driver fuelling the ongoing popula¬ 
tion decline, while habitat degradation, electro-fishing and 
fishing with poison also contribute to the species’ demise 
(Huang et al. 2008). The Vietnamese populations are cur¬ 
rently also threatened by habitat loss and alterations caused 
by intensive coal mining and illegal timber logging (Zieg¬ 
ler et al. 2008, M. van Schingen pers. obs.). The spe¬ 
cies’ small body size combined with its striking appear- 
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ance makes S. crocodilurus a desired target for poachers 
supplying the international pet trade (Le & Ziegler 2003, 
Huang et al. 2008). Therefore, the already heavily dimin¬ 
ished populations of S. crocodilurus will likely continue to 
decline in China as well as in Vietnam if no immediate pre¬ 
ventative conservation measures are initiated. While the 
species has not been assessed by the ‘lUCN Red List’ yet, 
it was assigned to Appendix II of the Convention on Inter¬ 
national Trade in Endangered Species of Wild Fauna and 
Flora (CITES) and classified as a Category-I species under 
the ‘Wild Animal Protection Law’ in China (Huang et al. 
2008, CITES 2013). This species has also been proposed to 
be included in the checklist of protected species in Viet¬ 
nam (Nguyen 2011). 

Correlative species distribution models (SDMs) have 
been used successfully to reveal potentially suitable habi¬ 
tats and investigate the effectiveness of protected areas (e.g., 
Araujo et al. 2004, Araujo et al. 2007, Hannah et al. 2007, 
Rodder & Schulte 2010, Rodder et al. 2010). The poikilo- 
thermic species’ strong dependence on environmental con¬ 
ditions (e.g., water, ambient temperature) (Zhao et al. 
1999, Ning et al. 2006, Wang et al. 2009) in combination 
with a preference for specific microhabitat characteristics 
renders S. crocodilurus an ideal taxon for performing spe¬ 
cies distribution modelling approaches. Thus, it is the aim 
of the present paper to predict the potential distribution of 
S. crocodilurus by applying SDMs and to identify potential¬ 
ly suitable habitats to guide further field exploration as a 
basis for improved protected area management planning. 

Methods 

We performed SDMs based on locality records and a set of 
environmental predictors that combine environmental var¬ 
iables and remote sensing data. We compiled a total of 20 
occurrence records, partly from our own field research in 
northern Vietnam and southern China, as well as from lit¬ 
erature (Huang et al. 2008). We computed a set of twelve 
environmental predictors based on temporal transforma¬ 
tions of remote sensing data, using the dismo and raster 
packages for Cran R (Hijmans & van Etten 2012, Hijmans 
et al. 2012, R Core Team 2012). A set of pre-processed re¬ 
mote sensing variables derived from MODIS sensors of 
two NASA satellites (spatial resolution = 30 arc sec; tempo¬ 
ral resolutions: MOD11A2 = 8-day averages; MCD43B4 = 
16-day averages [Mu et al. 2007, Scharlemann et al. 2008]) 
was obtained from the EDENext project (imagery pro¬ 
duced by the TALA Research Group, Oxford University us¬ 
ing methods described in Scharlemann et al. 2008). The 
raw remote sensing dataset comprised monthly averages of 
the enhanced vegetation index (EVI), the normalized veg¬ 
etation index (NDVI), and day- and nighttime land surface 
temperatures, collected between 2001 and 2005. The de¬ 
rived environmental predictors comprise variables describ¬ 
ing annual averages as well as seasonal variability (Tab. 1). 

We computed pairwise coefficients of determination 
based on Spearman rank correlations to assess co-lineari¬ 


ty. A subset of twelve variables was selected with R^ < 0.75, 
which was clipped to the spatial extent of the species’ geo¬ 
graphical range. We modelled the potential distribution of 
S. crocodilurus using the biomodi package v. 2.1.15 (Thuil- 
LER et al. 2013) for Cran R, applying the following algo¬ 
rithms: ‘Generalised Boosting Models’ (GBM), ‘Multivari¬ 
ate Adaptive Regression Splines’ (MARS), ‘Generalized Lin¬ 
ear Models’ (GLM), ‘Generalized Additive Models’ (GAM), 
‘Classification Tree Analyses’ (CTA), ‘Artificial Neuronal 
Networks’ (ANN), ‘Surface Range Envelopes’ (SRE), and 
‘Maxent’. The models were trained using a randomly select¬ 
ed subset of the species’ occurrence records (80%), while 
the remaining 20% were used to analyse model perform¬ 
ance with five iterations per algorithm, applying the ‘receiv¬ 
er operating characteristic curve’ (ROC) (Swets 1988), ‘Co¬ 
hen’s Kappa, and the ‘True Skill Statistic’ (TSS) (Allouche 
et al. 2006). We used 1,000 randomly created pseudo-ab¬ 
sence records within a circular buffer of 50 km, encom¬ 
passing each presence record for model building. Based on 
the SDMs, we computed an ensemble integrating all SDMs 
with ROC > 0.7 ranked according to their performance. 
The final ensemble was projected onto a rectangular area 
slightly larger than the area covered by the species records 
to highlight potentially suitable habitats in northern Viet¬ 
nam and southern China. As a presence/absence threshold, 
we selected the minimum score observed in the species’ 
records. Areas characterized by environmental conditions 
exceeding those available within the 50 km buffer enclosing 
all species records were excluded from projections as ex¬ 
trapolations beyond the training range of the ensemble, as 
these would likely increase the uncertainty factor. 

Potential habitat suitability for S. crocodilurus as well 
as the coverage with designated protected areas according 
to lUCN standards (categories I, II, IV, V, VI; lUCN 2013) 
were evaluated across the study area in order to reveal po¬ 
tentially suitable habitat yet unexplored for the occurrence 
of the species and to ease the future management planning 
of reserves. We obtained protected area shape files from the 
World Database of Protected Areas (lUGN, UNEP-WGMG 
2013). To characterize the realized and potential niche of 
the species, we extracted all environmental variables from 
the existing species records as well as from the available en¬ 
vironmental background within the 50 km buffer and com¬ 
puted density estimates using the sm package for Cran R. 

Results 

New population record 

As a result of our recent field research in northern Vietnam 
in 2013, another population of S. crocodilurus was discov¬ 
ered in the Dong Son - Ky Thuong Nature Reserve on the 
eastern side of Yen Tu Mountain in Quang Ninh Province. 
The new population is distant by about 40 km from the 
known subpopulations in Vietnam (Bac Giang and Quang 
Ninh provinces) and 380 km from the closest subpopula¬ 
tion in Ghina (Guangxi Province). We could also extend the 
known altitudinal range of the species in Vietnam (from 
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Table 1. Environmental variables and derived variables sets used for SDM development. 


Abbreviation 

Set variable 

Derived variable 

ED15078_biol0 

MODIS V4 Band 07 + 08 Synoptic Months: Day- + Night-time 
Land Surface Temperature 

Mean Temperature of Warmest Quarter 

ED15078_bioll 

MODIS V4 Band 07 + 08 Synoptic Months: Day- + Night-time 
Land Surface Temperature 

Mean Temperature of Coldest Quarter 

ED15078_bio2 

MODIS V4 Band 07 + 08 Synoptic Months: Day- + Night-time 
Land Surface Temperature 

Mean Diurnal Range (Mean of Monthly 
max-min Temp.) 

ED15078_bio3 

MODIS V4 Band 07 + 08 Synoptic Months: Day- + Night-time 
Land Surface Temperature 

Temperature Isothermality (BI02/BI07) (* 100) 

ED15078_bio4 

MODIS V4 Band 07 + 08 Synoptic Months: Day- + Night-time 
Land Surface Temperature 

Seasonality of Temperatures 

ED15078_bio5 

MODIS V4 Band 07 + 08 Synoptic Months: Day- + Night-time 
Land Surface Temperature 

Maximum Temperature Warmest Month 

ED15078_bio7 

MODIS V4 Band 07 + 08 Synoptic Months: Day- + Night-time 
Land Surface Temperature 

Temperature Annual Range (BI05-BI06) 

ED15I4_biol 

MODIS V4 Band 14 Synoptic Months: Normalised Difference 
Vegetation Index (NDVI) 

Annual Mean NDVI 

ED15I4_bio7 

MODIS V4 Band 14 Synoptic Months: Normalised Difference 
Vegetation Index (NDVI) 

Annual Range of NDVI 

ED15I5_biol 

MODIS V4 Band 15 Synoptic Months: Enhanced Vegetation 
Index (EVI) 

Annual Mean EVI 

ED15I5_bio5 

MODIS V4 Band 15 Synoptic Months: Enhanced Vegetation 
Index (EVI) 

Maximum Monthly EVI 

ED15I5_bio7 

MODIS V4 Band 15 Synoptic Months: Enhanced Vegetation 
Index (EVI) 

Annual Range of EVI 


400 to 800 m, see Le & Ziegler 2003, Ziegler et al. 2008) 
by discovering individuals occurring from 180 m a.s.l. in 
the Dong Son - Ky Thuong Nature Reserve to 850 m a.s.l. 
in the Yen Tu Nature Reserve, revealing a similar altitudi¬ 
nal range compared to Chinese populations (from 200 to 
1,500 m, see Huang et al., 2008, Zhao et al. 1999). 

Realized and potential niche 

With respect to all univariate comparisons, S. crocodilu¬ 
rus occupied an environmental niche nested in the over¬ 
all available niche space (Fig. 1). The comparisons revealed 
only slightly smaller spans of the realized niches com¬ 
pared to the available niches, but different density maxima 
(Fig. 1), as S. crocodilurus occupies areas with relatively low 
NDVI scores compared to the overall NDVI range, indi¬ 
cating its close dependence on intact vegetation. Further¬ 
more, S. crocodilurus inhabits areas with an annual tem¬ 
perature range that is relatively constant without extreme 
maxima. The overall relatively low annual temperature 
range is consistent with mountainous habitats. 

Potential distribution 

With ‘excellenf ROC values being obtained for the ensem¬ 
ble (ROC,^^,= 0.996, Kappa = 0.239, TSS = 0.957), the mod¬ 
el shows a strong ability to discriminate between suitable 


and unsuitable habitats. The variable EDi5i4_bio7 (Annu¬ 
al Range of NDVI) had the strongest effect on the model 
(52%) as measured by permutation importance, followed 
by EDi5078_bio7 (Temperature Annual Range) (50%), 
EDi5078_bio4 (Seasonality of Temperatures) (41%), 
EDi5078_bioio (Mean Temperature of Warmest Quar¬ 
ter) & EDi5078_bio3 (Temperature Isothermality) (38%), 
EDi5078_bio2 (Mean Diurnal Temperature Range), and 
EDi5i5_bio5 (Maximum Monthly EVI) (37%), whereas 
the remaining variables contributed on average less than 
35% to the final model. The ensemble revealed only small 
proportions of the study extent to provide suitable habitats 
(China: 1.12%; Vietnam: 4.29%; Fig. 2). Only a fraction of 
the selected study region was found to be covered by pro¬ 
tected areas in Vietnam (17.29%) while only 6.31% of the 
study region was found to be covered by reserves in China. 
Furthermore, only a fraction of habitat deemed suitable is 
presently located within designated protected areas (1.74% 
in China and 0.15% in Vietnam; Fig. 2). The model sug¬ 
gests potential additional occurrences of S. crocodilurus 
to exist, amongst others, in the Shiwandashan Nature Re¬ 
serve (SNR) in southern China. According to our mod¬ 
el, the SNR situated between the confirmed localities in 
southern China and northern Vietnam represents the larg¬ 
est contiguous area of potentially suitable habitat in China 
(Fig. 2). Small fragmented areas with high predicted prob¬ 
abilities for the occurrence of S. crocodilurus are scattered 
across northwestern Vietnam, but are presently not pro¬ 
tected. 
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Figure 1. Comparisons of density distributions of the realised bioclimatic space of Shinisaurus crocodilurus with the potential avail¬ 
able space along 12 environmental variables. Note that derived variables comprise relative scores specific to the study area. Therefore, 
only qualitative units are shown. 
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Discussion 

Density estimates across the 12 environmental gradients re¬ 
vealed different characteristics of the available vis-a-vis re¬ 
alised niche, indicating that S. crocodilurus is a habitat spe¬ 
cialist. Denoted already by Ning et al. (2006), the vegeta¬ 
tion index proved to be a determinant for the occurrence of 
S. crocodilurus. In accordance with previous studies (Zhao 
et al. 1999, Wang et al. 2009), temperature-related varia¬ 
bles revealed a strong contribution to the SDMs, as the spe¬ 
cies occupies habitats characterised by low temperatures as 
well as a low diurnal and annual temperature range. These 
microhabitat conditions are also characteristics of moun¬ 
tainous habitats. Our model revealed several spots covered 
by potentially suitable habitats to be situated in northwest¬ 
ern Vietnam. To date, the species has not been confirmed 
to occur in this area, and the Red River might serve as a 
geographical barrier, restricting the species’ distribution 
to northeastern Vietnam. However, SDMs are not able to 
identify geographical barriers and the accessibility of a po¬ 
tentially suitable habitat so that this hypothesis has to be 
verified by further field surveys. The small size of potential¬ 
ly suitable habitats combined with heavy fragmentation and 
poor coverage with designated protected areas underlines 
the urgent need for significant improvements of the existing 
reserve network to increase effectiveness for the conserva¬ 
tion of S. crocodilurus. Therefore, potentially suitable habi¬ 
tats with high detection probabilities should be surveyed 


for occurrences of the rare species. Due to its limited dis¬ 
persal capacity, its close dependence on water (Zhao et al. 
1999, Le & Ziegler 2003, Zollweg & Kuhne 2003, Zheng 
& Zhang 2004, Ning et al. 2006, Huang et al. 2008), and 
a rather sedentary lifestyle, we expect locations situated in 
close proximity or between confirmed populations, such as 
the SNR in southern China or the Khe Ro Sector within 
Tay Yen Tu Nature Reserve in Vietnam, to be most promis¬ 
ing. Such areas might represent important stepping stones 
for the species. The few existing reserves presently hold¬ 
ing populations of S. crocodilurus (Tay Yen Tu, Yen Tu, and 
Dong Son - Ky Thuong nature reserves in northern Viet¬ 
nam and the Luokeng, Daguishan, and Linzhouding nature 
reserves in southern China) need to be subjected to signifi¬ 
cant law enforcement to reduce the collection of individuals 
to a minimum and prevent electro-fishing and fishing with 
poison. Moreover, these protected areas should be consid¬ 
ered for status elevation to prevent further habitat loss and 
fragmentation. The alarmingly rapid population declines 
observed recently throughout the species’ distributional 
range highlights the urgency of an assessment for the lUCN 
Red List, which is currently undertaken by the lEBR and 
the Cologne Zoo, as well as the need for a zero quota on 
the commercial trade of wild-caught specimens. Not only 
the status, but also the size of the remaining populations 
should be analysed or re-analysed contemporarily to clarify 
whether scientifically coordinated population restoration is 
required besides improved habitat protection measures. 



Figure 2. Occurrence records of Shinisaurus crocodilurus are displayed as black circles, with potential habitat suitability ranging from 
low (yeUow) to high (red), coverage with designated reserves (stippled polygons), and the course of Red River (blue). For dark grey 
areas, no predictions could be made, as environmental conditions exceed the training range of the SDM. Only vague locality informa¬ 
tion is displayed in order to protect remnant populations. 
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